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^ ' We consider the dynamics of systems with arbitrary friction and diffusion. Tliese include, as a spe- 

O \ cial case, systems for which friction and diffusion are connected by Einstein fluctuation-dissipation 

, relation, e.g. Brownian motion. We study the limit where friction effects dominate the inertia, i.e. 

where the mass goes to zero (Smoluchowski-Kramers limit). Using the Ito stochastic integral con- 
vention, we show that the limiting effective Langevin equations has different drift fields depending 
on the relation between friction and diffusion. Alternatively, our results can be cast as different 
interpretations of stochastic integration in the limiting equation, which can be parametrized by 
(~| ' a € R. Interestingly, in addition to the classical Ito (a = 0), Stratonovich [a — 0.5) and anti-Ito 

O ^. (q — 1) integrals, we show that position-dependent a = a{x), and even stochastic integrals with 

I i Q ^ [0, I] arise. Our findings are supported by numerical simulations. 

I. INTRODUCTION 

^ , Most physical, chemical, biological and economic phenomenapresent an intrinsic degree of randomness. These are 

^ • typically modelled by stochastic differential equations (SDEs) [1| . SDEs were introduced at the beginning of the XX 
' century to describe Brownian motion by adding a random driving function to an ordinary differential equation (ODE); 
, since then, SDEs have come into widespread use in, e.g., physics, biology, and economics. However, SDEs involving 
■ multiplicative noise terms can be integrated according to various definitions leading to different solutions 0], e.g. 
, the Ito integral and the Stratonovitch integral. From the modeling point of view it is, therefore, key to know what 
04 ' definition to use in any given situation Q . From the mathematical point of view the simplest approach is to write all 
equations according to, e.g., the Ito definition of the stochastic integral; the different interpretations mentioned above 
then reappear as additional drift terms, which are often referred to in the literature as "spurious drifts." 

In order to understand the origin of the difficulty, we consider a simple case. The SDE dxt = a{xt)dWt, where Wt 

is a Wiener process, can be solved by integration, i.e. Xt = xo + / a-{xs)dWs, where we define / cr{xs) o" dWg = 

Jq Jq 

u ■ ^ 

C3 ■ lim a{xt»)AWt„ with AWt„ = Wt„.^ - Wt^, <„ = f and = ^f^t and, typically, a G [0,1]. Since Wt is 



a function of unbounded variation, differently from ordinary Riemann-Stieltjes integrals, the limit of these partial 
summations generally leads to different values of the integral depending on the choice of a. In particular, a = leads 
to the Ito integral, a = 0.5 leads to the Stratonovitch integral, and a ~1 leads to the anti-Ito (or isothermal) integral. 
This is the reason why it is necessary to give both a SDE and the respective a with which to solve it in order to have 
a fully determined model Q. 

Various preferences regarding the appropriate choice of a have emerged in the numerous fields where SDEs have 
been a ppliedQ. For example, the martingale property, i.e. the specific feature of the Ito integral of "not looking into 
the future," meaning that, when the integral is approximated by a summation, the leftmost point of each interval is 
used, is the main reason of its popularity in economics [l| and biology 0, Q. In general, the Stratonovitch integral 
emerges naturally when the Wiener process is replaced by a sequence of approximating deterministic processes and 
has the advantage of leading to ordinary chain rule formulas under a change of variable Q. However, the fact that 
Stratonovitch integrals are not martingales gives the Ito integral an important computational advantage Finally, 
the anti-Ito integral has been shown to be the most appropriate to describe ph ysic al phenomena that are in e quilibrium 
with a heat-bath for which Einstein fluctuation-dissipation relation holds |9l4l3l|. In particular, equations that satisfy 
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FIG. 1: (a) For a Brownian particle 7(2:) (dark line) and a-{x) (grey line) are related by the Einstein fluctuation-dissipation 
relation [Eq. Q]; in this case, 7(2:) = (1 -fa;/100) (b) The solution of the Newton equations [Eq. ([l}] for m — >■ (dashed lines) 
converge to the solution of the SK approximation [Eq. ((2}] for a = 1, i.e. anti-Ito integral, (black solid line); also the solution 
for a = 0, i.e. Ito integral, (grey solid line) is given for comparison. All solutions are obtained for the same Wiener process. 
The inset is a blow-up of the flnal part of the trajectories (dashed square). 



the fluctuation-dissipation relation occur in molecular dynamics. Here, the hmiting equation is constrained to be the 
anti-Ito type to correct the invariant distribution to model the Gibbs distribution. 

A classical example of a phenomenon in equilibrium with a heat bath is the motion of a mesoscopic particle of mass 
m immersed in a fluid, i.e. Brownian motion. If we assume that the particle moves in one dimension under the action 
of a continuous force F(x), its position being a;™ G R at all times t > in a finite interval, the corresponding Newton 
equation is: 



dec" 



mdvY" = F{x];')-j{xY')vY'dt + a{xT)dWt 



(1) 



with the initial conditions = v™ and x™' — x'"- The friction coefficient 7(2;) > and the intensity (standard 
deviation) of the noise <y{x) > are, in general, position-dependent; we also assume that F, 7, and a are differentiable 
functions of x smooth enough so that the process {xt,Vt) exists for all i on a finite interval. It is well known that, 
since the derivative of Xt, i.e. = w™, exists, the stochastic integral in Eq. ([1]) is equivalent under all interpretations 

The limit of Eq. ([l} as jti — >■ has to be studied with care, requiring a nontrivial computation. In general, similar 
limits involve additional drift terms, i.e. "spurious drifts." A statement of our results in terms of different definitions 
of stochastic integral is also possible and, in some cases, straightforward. By naively setting m = in Ec^. ([T]), we 
obtain a SDE for — Xt'. 



dxt 



dt 



cr{,xt) 



dWu 



deterministic stochastic 



(2) 
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with initial condition xq = x^- Eq. ([2|) is called the Smoluchowski-Kramcrs (SK) approximation to Eq. ([T]). Differently 
from the solution of Eq. ([T]), the solution of Eq. ^ depends on the interpretation of the stochastic term, i.e. on the 
choice of a. In the following we identify the value of a that introduces the correct additional drift; we emphasize that 
this is just another way of expressing additional drift terms within the Ito formalism. 

We can gain some insight into this zero-mass limit procedure considering numerical solutions of Eq. ([T]) for various 
decreasing values of m, but for the same Wiener process [Fig. [Tj. For a Brownian particle the Einstein fluctuation- 
dissipation relation holds: 

j{x) (X a{x)'^. (3) 

In Fig. [ija) 'y{x) (dark line) and a{x) (grey line) are presented. The dashed lines in Fig. \V[h) represent some solutions 
of Eq. ([1]) for decreasing m: they become rougher and rougher as the m decreases. They converge towards the anti-Ito 
(a = 1) solution of Eq. (0_(black solid line); this is in agreement with the recent experimental demonstration [T3| 
and mathematical proof that for a Brownian particle the most natural interpretation is the anti-Ito integral. We 
remark that the Ito {a = 0) solution of Eq. ([2]) (grey solid line) presents clear deviations from the correct one, as can 
be clearly seen in the inset of Fig. [H^b). 

In this article we study the zero-mass limiting behavior of a larger class of equations that have the form of Eq. ([1]), 
but for which 7(0;) and cr(x) are allowed to var y in dependently from each other. This can be the case, e.g., in the 
description of the evolution of complex systems [la] . We study this general class of equations from the point of view 
of the convergence of the infinitesimal operators of the corresponding diffusion processes , i.e. backward Kolmogorov 
equations, thus following the well-known methods from homogenization theory [16l4l8l |. The main result of this paper 
identifies for given ^{x) and cr(x) the drift term and hence the corresponding a [Eq. (jl9p ]: we find that in general it can 
be a function of x. We remark that the spurious drifts are defined assuming the Ito stochastic calculus convention, 
while the values of a are defined with reference to SDE ([J); we introduce the notation o"(^)c?M^( to indicate the 
presence of such extra drift. Interestingly, we find that when a generalized fluctuation-dissipation relation holds, i.e. 

-f{x)(xa{x)\ (4) 

a is only a function of the exponent A and independent from x [Eq. (j22p ]. In particular, for A = wc retrieve the 
Ito interpretation and for A = 2 the anti-Ito interpretation, while the Stratonovich interpretation is only retrieved 
asymtotically for A — )• 00. Interestingly, also values of a ^ [0,1] occur for A € (0,2). Although in this article we 
always consider the variable x to be one dimensional, the general multi-dimensional case can be studied using similar 
methods [l7j . 

In section [Til we give further details on the ansatz we will use. In section IIIIl we introduce a relation [Eq. ([5]) 
and Eq. ©J that permits one to express any stochastic integral as an Ito integral with a modified drift term. In 
section IIVI we analyze an asymptotic expansion of the solution of the Kolmogorov equation in order to derive the 
effective SK-approximation [Eq. ([T5)) ]. In section |Vl we derive a formula for the value of a{x) as a function of j{x) 
and (t(x) [Eq. (fTO)) ] and we study various important special cases: the case of constant jix) [? IV A| : the case of a 
Brownian particle for which Eq. ([3]) holds [S IV B| : the cases for which Eq. holds, which are all the cases where a is 
position-independent [§|VC]; and the singular case for which 7(0;) (x a{x), which leads to a SDE without ambiguity, 
but with an additional drift nevertheless [ WD| . 



II. THE ANSATZ 



Our ansatz is to study the zero-mass limit of a diffusion process defined by Eq. ([1} with arbitrary 7(0;) and 
(t{x), analyzing the behavior of its infinitesimal generator, i.e. the Kolmogorov equation. We perform an asymptotic 
analysis of the Kolmogorov equation, expanding its solution in powers of a small parameter, i.e. ^/m. The result is a 
drift term in the effective SDE, which can be translated into the correct interpretation of the stochastic integral. In 
order to clearly present the results in a technicality-free style, the approach is not fully rigorous. Similar convergence 
results have been studied [l9j . 

The problem of taking the mass to zero, justifying the limiting equation, has been addressed by various authors 
at different levels of generality and mathematical rigor, starting with M. Smoluchowski [23| and H. Kramers [llj. E. 
Nelson studied the F — case in which 7 and a were constant and proved that the solution of Eq. ([T]) converges to 
the solution of Eq. ([5]), which in this case was unambiguous since a was constant [22|. The case with constant a but 
including an external force was also treated (see reference [l^ and references therein) with a similar result but by 
entirely different methods. The case where a depended on the position, but 7 was constant, was first studied by M. 
Freidlin in [23j . showing that the limiting equation should be interpreted with a = and that, in the presence of the 
colored noise, in some cases the limiting a equals ^. Subsequently, R. Kupferman et al. showed that all values of a 
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between and ^ could be obtained by taking the correlation of the noise and the mass of the particle to zero in an 
appropriate way [13] . 

III. STOCHASTIC INTEGRALS AS MODIFIED ITO INTEGRALS 

A stochastic integral with a given a{x) can be expressed as an Ito integral, i.e. a = 0, with an additional noise- 
induced drift, i.e. a "spurious drift." To justify this claim, we consider an arbitrary process 

dxt - b{xt)dt + a{xt) o"(=") dWt, 

where the SDE is defined by a{x) and the solution xt, is real valued and one dimensional. The integrated equation is 

N 



Xt ^ Xo + 



/ b{xs)ds+ lim V cr(a;t=)AW^t„, 



with = By expanding (j{xt), we see that this corresponds to 

/'* /'* da(x ) ^ 

xt^x^ + I b{xs)ds+ / a{xs)a{xs) " ds + lim (T{xt^)AWt„, 
Jo Jo dxs 

with tn — ^t. And this can be interpreted as the Ito (a = 0) equation 

dxt = b{xt)dt + a{xt)a{xt)—j-^ + a{xt)dWt, (5) 

dxt 

where we omit the o'^ for all further Ito integrals. In particular, Eq. ([2]) interpreted with any a{x) corresponds to the 
Ito equation: 



dxt 

with initial condition xq = x"- 



F{xt) , , . cr(xt) d f cr{xt) 
a{xt) 



j{xt) "fixt) dxt \"fixt) 



dt+^dWu (6) 



IV. SMOLUCHOWSKI-KRAMERS APPROXIMATION 

In order to simplify further analysis, we substitute u™ — y^w"' in Eq. ([T]) obtaining the following two-dimensional 
SDE: 

C dxj" = -^wj" dt 



— * ut 



dt 4- dWf 



with initial conditions = x™ ^-^d u™ = y/miy"^. 

To determine a we use a multiscale analysis of the backward Kolmogorov equation of the SDE (O. Let 
g{x' ,u' ,t'\x,u,t) be the probability density of the distribution of the position and (rescaled) velocity {x',u') of the 
particle at time t' given their values (x, u) at a time t < t' . Then the backward Kolmogorov equation for the SDE ([7]) 
is 

dg{x' ,u' ,t'\x,u,t) a{x)^ d^g{x' ,u' ,t'\x,u,t) u dg{x' ,u' ,t'\x,u,t) 

- ' — (8) 



dt 2m dv? ^Jm dx 

F{x) dg{x' , u' , t'\x, u, t) "f{x)u dg{x' , u', t'\x, u, t) 
du m du 



Since the equation involves derivatives with respect to the a:, u, and t variables we write g{x' , u' , t'\x, u, t) = g(x, u, t), 
to shorten notation. Eq. can be rewritten as 

— = —Li + -^Li g. 

Ot \m Jm 
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with 

d d 

Notice that the operator Li is the generator for an Ornstcin-Uhlcnbeck (OU) process with coefficients dependent on 
X. We denote this process as u, and write the stochastic differential equation 

dut ~ —^{x)ut dt + a{x) dWt- (9) 

The invariant density for ut is 

g*(^)^C(x)exp(^ ~j|^jf ^ , (10) 

where C{x) is a normahzing constant. We postulate that the solution of the Kolmogorov equation has an asymptotic 
expansion 5 = (?o + V^gi + rn92 + ■■■ IWm. We match powers of m to obtain the following equations, 

Ligo = 0, (11) 
iiffi = -^2.90, (12) 

^ = ^1.92 + ^251, (13) 
where Li,L2 are differential operators. Solving Eq. (|lip results in 

go{x,u,t) ^ Ci{x,t) du + C2ix,t). 



Since the first term is not integrable in u, Ci must be zero and thus go ~ go{x, t), independent of u. By the Fredholm 
alternative, the solvability condition for Eq. ([T2|) is given as 



/oo 
g*L2go du = 0, 
-00 



for all 5* such that Llg* = [25|, l26| . Here L* is the adjoint of Li: 



The relevant (integrable) solution is a mean-zero Gaussian given in Eq. (jlOp . which satisfies g*{u) = g*{—u), thus 

[ 9* L2go du = f ug*{u)du^O. 
Next the solvability condition for Eq. ([T^ is 

p I -L2L^'L2go + ^^g*du = 0. (14) 



First we set 



V = L^^L2go, 



which by the previous solvability condition is well defined. Thus LiV = ^230: or 

aix)^ d^V , . dV _ dgo 



2 du'^ du dx 

Notice that the function 

y ^ u dgo 

j{x) dx 



'y{x)u—— = u——. (15) 
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is a particular solution of Eq. p5|) . From Eq. p4)) we must have 



a^^o , d-f{x)dgo F{x) dgQ\ dgo 

du^- — , (16) 



~7(x) dx^ 7(2;)^ ax 7(2;) dx J 5t ' 

for any satisfying L\g* = 0, in particular for the invariant density of the OU process ut- After Gaussian integration 
over u Eq. (|16p becomes 

a{xf d^go , f F{x) a{x)^ d'y{x)\ dgo _ dgo ^^^^ 



2^(x)^ dx^ \ 7(2;) 2^{x)^ dx J dx dt 
This gives the SK approximation to the backward Kolmogorov equation. The corresponding (Ito) SDE is 

= - + ^ dW.. (18) 

\"t{xt) 2'y{xtr d,x J '-f(xt) 

Since we derived this equation from the convergence of the infinitesimal operators, rather than directly studying the 
limit of the SDE ([7]), the convergence is in law. 

V. AN EQUATION FOR a{x) 

We derive an equation for 01(2:) depending on the friction j{x) and diffusion (t{x) comparing the backward Kol- 
mogorov equation of the SDE [Eq. ^] and Eq. and solving for a{x): 

a(x) = I'i^M^) .19^ 

^ ' 2(7'(x)a(x)-7(x)a'(cc))' ^ ' 

where 7'(x) = ^^^^ and cr'(x) = This equation shows that in general, a varies with position and can even 

take values outside the interval [0,1]. Interestingly, a never takes the value i, i.e. we never obtain a Stratonovich 
correction. 

A. 7(2;) = 70: Constant friction 

The case in which ■y(x) = 70, while <j{x) is allowed to vary [Fig, dja)], has been often object of mathematical studies. 
For example, Freidlin |23| and later Pavliotis and Stuart [27[ proved that the limiting equation has a stochastic term 
with a = 0; this result is rederived here. Physically, an example of this system is in the framework of the Maxey- Riley 
model of inertial particles in a Gaussian field [2^ with correlation time assumed to be very short. In Fig. [2jb), we 
show how the numerical solutions for m — >■ converge towards the Ito (a = 0) solution of Eq. 

B. 7(2:) oc a(x)'^: Brownian motion 

The particular case when Einstein fluctuation-dissipation relation is satisfied in its standard form [Eq. ([3])] is 
particularly important because it describes the diffusion of Brownian particles. If D{x) denotes the hydrodynamic 
diffusion coefficient and ksT the thermal energy, then 

7(^) = ^ (20) 
D{x) 

and 

^/D{x) 

This case was studied experimentally in [T2|. [l3|. showing that the correct value of a for m — > is a = 1. It was 



subsequently proven by a mathematical argument jl4| that in this case the processes a;™ converge to the solution of 
the limiting equation with a = 1 in the sense; this result is rederived here, but only in a weaker sense. 
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FIG. 2: (a) a-{x) (grey line) and 7(2;) = o-{x)° = constant (dark line), (b) The solutions of the Newton equations [Eq. fl}] for 
m ^ (dashed lines) converge to the solution of the SK approximation [Eq. ((2|] for a = 0, i.e. Ito integral, (grey solid line); 
also the solution for a = 1, i.e. anti-Ito integral, (black solid line) is given for comparison. All solutions are obtained for the 
same Wiener process. The inset is a blow-up of the final part of the trajectories (dashed square). 

C. ^{x) oc a{x)''^: Constant a 



All the cases for which a{x) = a, can be obtained equating the right-hand side of Eq. (fT9|) to a constant, different 
from i. After a simple calculation, we obtain 7(0;) = ca{x)^, where c is a constant, i.e. Eq. Qj. It follows that 



A 



2(A-1)- 



(22) 



The value of a depends on the exponent A as shown in Fig. [3l This result includes as particular cases 7(.t) = 7, for 
which a = [§ IV A| . and j{x) = C(t(x)^, for which a = 1 [S IV B| . However, we remark that the value a = ^ is only 
obtained asymptotically for A — >■ cxj. 

Interestingly, values of a outside the interval [0, 1] can be achieved for certain friction-diffusion relations. For 
example, the relation 7(2;) = a{xY^'^ gives a = 2 by the formula [Eq. ([22]) ]. Figure (j4]) gives insight to the zero mass 
limit. Different constructions of the stochastic integral are given for Ito (grey solid line), anti-Ito (black solid line), 
and for a = 2 (dark grey solid line), for the same Wiener process. 



D. 7(2;) cx a{x): A singular case 



When 7(x) cx a{x) [Fig. [5l^a)], the stochastic term in Eq. ([2]) gets multiplied by a constant factor, i.e. 



constant, and thus there is no ambiguity in its solution. However, the zero-mass limit of Eq. ^ does not converge 
to this solution. This can be seen by setting j{x) ~ ca{x) and using Eq. (|17l) directly; the limiting equation is 



dgo 
dt 



1 3^50 , ( F{x) a'{x) 



2c^ dx"^ 



c<t{x) 2c^a{x) 



dgo 
dx 



(23) 
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X 

FIG. 3: a as a function of A for the case when 7(2;) oc o"(a;)'*' [S IV Cj . For A — > l,a diverges asymptotically (dashed line) leading 
to the singular case discussed in § IV Dl The Ito integral (a — 0) is obtained for A = (square) and the anti-Ito (a — 1) for 
A = 2 (diamond); the Strasonovich integral (a — 0.5) is only obtained asymptotically (dotted line) for A — 00. 

This gives the Smoluchowski-Kramers limiting SDE as 

dt+'-dm. (24) 

\ca{xt) 2c^a{xt) J c 

Here we see that there is a correction to the drift term. In Fig. El^b), we show how the numerical solutions for m — > 
(dashed lines) converge towards the solution of the SK approximation (black solid line), while the solution Eq. ([2]) 
without the correction to the drift clearly diverges (grey solid line). 

VI. CONCLUSION AND FUTURE RESEARCH 

We have performed a systematic study of the Smoluchowski-Kramers limit for a class of SDEs with arbitrary 
friction and diffusion. We have identified the Ito form of the limiting Langevin equation and discussed its equivalent 
interpretation in terms of other definitions of stochastic integrals. We introduced the notation o"(^) dWt, that is 
interpreted as an Ito integral with an additional well-defined drift term. However, we have not given a rigorous 
mathematical construction of this integral, which is nevertheless of interest. In a future work, we will study the 
analogous problem for some physically relevant systems driven by colored noise. 
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FIG. 4: (a) o"(a;) (grey line) and 7(2;) — a(x)'^^^ (dark line), (b) The solutions of the Newton equations [Eq. ([T])] for m ^ 
(dashed lines) converge to the solution of the SK approximation [Eq. ([2])] for a — 2, (dark grey solid line); also the solution for 
a — 0, i.e. Ito integral (grey solid line), and a — 1, i.e. anti-Ito integral, (black solid line) is given for comparison. All solutions 
are obtained for the same Wiener process. The inset is a blow-up of the final part of the trajectories (dashed square) 
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